CONVERGENCE OF SINGULAR LIMITS 
FOR MULTI-D SEMILINEAR HYPERBOLIC SYSTEMS 
TO PARABOLIC SYSTEMS 

DONATELLA DONATELLI AND PIERANGELO MARCATI 

Abstract. In this paper we investigate the zero-relaxation limit of the fol- 
lowing multi-D semilinear hyperbolic system in pscudodiffcrcntial form: 

W t (x,t) + -A(x,D)W(x,t) = —B{x,W(x,t)) + -D{W(x,t)) + E(W(x,t)). 
e e z e 

We analyse the singular convergence, as e J, 0, in the case which leads to a limit 

system of parabolic type. The analysis is carried out by using the following 

steps: 

(i) We single out algebraic "structure conditions" on the full system, moti- 
vated by formal asymptotics, by some examples of discrete velocity mod- 
els in kinetic theories. 

(ii) We deduce "energy estimates " , uniformly in e, by assuming the existence 
of a symmetrizer having the so called block structure and by assuming 
"dissipativity conditions " on B. 

(iii) We perform the convergence analysis by using generalizations of Com- 
pensated Compactness due to Tartar and Gerard. 

Finally we include examples which show how to use our theory to approx- 
imate prescribed general quasilinear parabolic systems, satisfying Petrowski 
parabolicity condition, or general reaction diffusion systems. 

Key words and phrases: Hyperbolic systems, parabolic systems, pseu- 
dodifferential operators, relaxation theory. 



1. Introduction 

In this paper we study the following semilinear multidimensional hyperbolic sys- 
tem with a small parameter e > 

(1.1) W t (x,t) + ^A(x,D)W(x,t) = -^B(x,W(x,t)) + ^D(W(x,t)) + E(W(x,t)), 

where W = W(x,t) belongs to R N , x E R d , t > 0, A{x,D) is a first order 
pseudodifferential operator. The system (1.1) includes also the case of hyperbolic 
differential operators of the form 

1 d 1 1 

(1.2) W t (x,t) + - A J (x)d J W(x,t) = — B(x, W(x,t)) + -D(W(x,t)) + E(W(x,t)), 

3=1 

where Aj(x), j — 1, . . . ,d are N x N matrices for any x C M. d . Our aim is to 
describe the limiting behaviour of the system (1.1) as e goes to zero. We look 
for structure condition to ensure that (1.1) approximate a second order parabolic 
system. Our interest in this problem is motivated also by a very strong similarity 
with the limiting structure appearing in the investigation of the hydrodynamic 
limit for the Boltzmann equation, in particular in the discrete velocity case. The 
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Boltzmann equation describes the evolution of the density f(x,^,t) of particles 
which are at time t in position x with velocity £ and has the following form 

vft + &-f = - e Q{f,f) 

where v is the Mach number and e the Knudsen number. By averaging f(%,£,t) 
in £ and by using higher order momenta we can define a hierarchy of macroscopic 
quantities. The investigation of the hydrodynamic limit regards the behaviour 
of those quantities (actually combined with the closure problem) as the Knudsen 
number goes to zero. In the case where the Mach number is of the same order 
of the Knudsen number our limit can be described by the Navier Stokes equation, 
otherwise when the Mach number is fixed and the Knudsen number tends to zero we 
end up with the Euler equation. Those limits allow us to understand the differences 
between relaxation limits of hyperbolic type to parabolic with respect to those one 
of hyperbolic to hyperbolic type @, @, @> Il> II > & & In particular if 
we deal with a discrete velocity models, the equivalent " Boltzmann equation" is 
a semilinear hyperbolic system. The simplest example is given by the Carleman's 
equation 

I 1 2 2 

fit + - fix = ^ (/ 2 - fi ) 

fit f2x — -o(/l _ /!): 

where we take £ G { — 1,1} and f\ — f(x,l,t), fi — f{x,—l,t). By rescaling the 
variable we get that p = f\ + fa as e [ satisfies the nonlinear diffusion equation 

Pt = ^OgpW 

This asymptotic problem was first investigated by Kurtz [^lj and McKean p6[ . 
Therefore the nonlinear diffusion problem, obtained as the limit of the Cattaneo 
hyperbolic nonlinear heat conduction equation, was proved by Marcati, Milani and 
Secchi |32| . The paper of Marcati and Milani |HJ concernes the pourous media flow 
as the limit of the Euler equation in 1-D, later generalized by Marcati and Rubino 
|35[ to the multi-D case. Relaxation phenomena of the same nature appear in the 
zero relaxation limits for the Euler-Poisson model for the semiconductor devices 
and it was investigated by Marcati and Natalini |33| , Q in the 1-D case and by 
Lattanzio and Marcati J2^] in the multi-D case. More recently Lions and Toscani 
(27) investigated a discrete velocity model leading to the pourous media flow. All 
of these papers, with exception of 0), [[36), make use of the techniques of compen- 
sated compactness. Similar ideas have been applied by Marcati and Rubino |}5| to 
show the general theory for 2 x 2 systems in the 1-D case and to propose a gen- 
eral framework that we are going to investigate here in the semilinear system case. 
Models of BGK type approximation have been successfully studied in this frame- 
work by Bouchut, Guarguaglini and Natalini |l| and Lattanzio, Natalini (2^] in the 
case of 1-D systems. The framework of |35| was also investigated in the quasilinear 
case by Lattanzio and Yon g for H s - smooth solutions. Recently, with a similar 
approach, Junk and Yong ]15[| derived the incompressible Navier Stokes equations 
form the BGK model. Preliminary results concerning semilinear systems have been 
obtained in , S , in particular in [S we considered a 1- D semilinear system with 
variable coefficients. Already in that case the classical compensated compactness is 
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not sufficient and it is necessary to use a generalization of this theory due to Tartar 
|46| and P.Gerard [jfl). From the technical point of view the problem, here, has 
additional complications , since we deal with a multi- D pseudodifferential system 
and since the symmetrizers are pseudodifferential operators. 

The plan of the paper is as follows. In Section 2 we give the basic definitions, the 
basic notations and we recall some mathematical tools needed in this paper. In 
Section 3 we describe the structure condition on our system and on its symmetriz- 
ers and we will describe the formal framework of the limiting process. The Section 
4 is devoted to give rigorous proof of the previous formal analysis, the assumptions 
of the previous section will allow us to obtain the a priori estimates of energy type 
uniformly in e. By using only the informations coming from the energy estimates 
and by means of the previous mentioned compactness framework we will be able 
to obtain our relaxation results. The structure conditions used in the limit process 
allow us al so t o satisfy the parabolicity condition for the limit system in the case 



of system (L2) . In Section 5 we show how this theory can be greatly simplified in 
the case of constant coefficients. Finally, in Section 6 we show how it is possible to 
use the theory developed in the previous sections in order to approximate a given 
parabolic system by means of a larger hyperbolic system with simpler nonlinear 
structure. We are able to approximate the following nonlinear parabolic systems in 
divergence form 

U t + J2 °i ( F *W - E BtiMOiUj = G(U) 

where x £ R d , t £ R+, U — U(x,t) £ R fc . With the previous techniques we 
approximate also reaction- diffusion systems of the form 

d 

u t=J2 A jA*)djd k U + f(U) 

where x £ R d , t £ R + , U = U(x, t) £ R fc . In this latter case we use two different 
approach. The former is based on the theory of pseudodifferential operators while 
the latter is based on symmetric differential operators and can be more usefull for 
numerical computations. Let us remark that this latter approximation extends 
in multi-D those proposed by Jin and Liu jlj] and Lattanzio and Natalini | j23| 
(example 5.3). Our theory does not require the use of L°° estimates (for instance 
via invariant domains like Serre p3[ ). 

2. Preliminary notions 

We start this section introducing the main notations and definitions used in this 
article. In particular we recall some basic facts and notations concerning the theory 
of pseudodifferential operators and we also recall our main compactness tools used 
in the strong convergence analysis. Namely 

(a) (•, •) denotes the scalar product in R 9 , (q — 1,2, ...) and | ■ | the usual norm of 

R ? (9 = 1,2,4 

(b) Ai m xn denotes the linear space of m x n matrices, 

(c) (•, -) 2 denotes the scalar product in L 2 (R d ) and || • || the norm in L 2 (R d ), 

(d) D(R d xR + ) denotes the space of test function C 00 (K d xR + ), V'(R d xR+) the 
Schwarz space of distributions and (•, ■) the duality bracket in V(R d x R+), 
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(e) H is a separable Hilbert space, C(H) the space of bounded operators, IC(H) 
the space of compact operators, 

(f) we denote by Hf (Q,H) the classical local Sobolev space of order s, i.e. 
u G fff oc (0, H) V v G C °°, (05) g L 2 (R", (1 + |£| 2 ) s d£). 

(g) we denote by F^Z 1 , Z 11 ) the derivative respect to the variable Z 11 oiF(Z I \Z If ). 

We shall also make use of the notion of parabolicity for systems of equations in 
various way (see Taylor |l9| volume III, [Q, Eidel'man Kreiss and Lorenz pc|). 
Let us consider the following system 

(2.1) ut = ^2A j ' k [t,x,D' x u)d j d k u + B{t,x,D' x u), 

j.k 

where u G W, A j ' k (t,x,D' x u) G M pxp , B(t,x,D' x u) G W and D' x is a differential 
operator of order not greater than two. The system is said strongly parabolic if 
there exists cq > 0, such that for all £ G R d one has 

J2A j ' k (t,x,D' x u)^ k >c |£| 2 /. 

Namely, if we denote L(t, x,D' x u, £) = — A J,fc (t, x, D' x u)£j£k this condition is 

equivalent to say L + L T is a negative definite matrix. Unfortunately, this con- 
dition is often difficult to be verified, then we recall a more general notion often 
referred as Petrowski par abol icity (see Taylor 49 volume III, J48fl). 
We say that the system (2.1) is parabolic if, denoted by Xk(t,x, D' x u,^) the eigen- 
values of the matrix L(t,x, D' x u,£), one has there exists ao > such that, for all 
£ G R d , 

ReX k {t,x,D' x u,0 < -aolCl 2 . 

The latter notion of parabolicity is equivalent to ask the existence of a symmetric 
matrix Po(t, x, D' x u, £), positive definite (i.e. Pq > cl > 0), homogeneous of degree 
in £, such that 

-(P L + L*P ) >c|e| 2 />0. 

Let us recall the basic notations concerning pseudodifferential operators to be used 
later on, we refer to |48| for details. Assuming p, S G [0, 1], me R, we denote S™ s 
the set of C°° symbols satisfying 

\Df>D%p(x,t)\ <C a<f3 (0 m - plal+m 

for all a, 0, where (£) = (1 + ICI 2 ) 1 ^ 2 - In such case we say that the associated 
operator defined by 

P(x,D)f(x) = / p(x,$f(0e ix tdZ := OP(p{x,Q) 



(where/(£) = (2tt) n J f(x)e ix ^dx denotes the Fourier transform of the function 
/) belongs to OPS™ s . If there are smooth symbols p m -j(x. £), homogeneous in £ 
of degree m-j for |£| > 1, i.e. p m ^(i,r() = r m ~ p m -j (x, £) for r, |£| > 1, and if 

j>0 
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in the sense that 

N 
j>0 

for all N, then we say p(x,£) E S m . We define also 

5,-00 = n 

m>0 

The following properties will also be used here. 

Theorem 2.1. (@, Theorem 0.5 A) If P(x,D) E OP5° 5 and < <5 < p < 1 

P(x,i?) : £ 2 (R") — > L 2 (R n ). 



Theorem 2.2. ((47), Chapter II, Lemma 6.2 Ex.8.1, Chapter III, Proposition 3.1) 
If P(x,D) E OPSp g , S < p and ifUep(x,£) > c > i/ien there exists R(x,D) E 
OPS® s such that R(x, D) > r]I > 0, KeP(a;, D) = ±(P + P*) 7 

3?eP(ir, £>) - P(x, D) E OPS~°°. 



Proposition 2.3. (Commutator estimate) ( p8| , Proposition 3.6 B) 
Given P(x, D) E OPS® we Ziaue 

l|[P, /]«||w.«> < c||/||iip||it||w"- /or0<<j<l, 
m particular if a — and p = 2 

ll[P, /]"IU 2 (R") < c||/||« P ||u|| H -i. 

Finally we state here our main tool, due to Tartar and Gerard ( p6| , Ql)), to 
study the convergence of quadratic forms with variable coefficients. The classical 
results concerning the use of Compensated Compactness in the theory of hyperbolic 
systems are reported in the books of Dafermos || and Serre |42|| . Let H, i?" denote 
separable Hilbert spaces, ft E K", an open set, m E N. We have the following 
theorem taken from Gerard jll] . 

Theorem 2.4. (Compensated Compactness) 

Let P E OPS m with principal symbol p(x,£) and {uk} be a bounded sequence of 
Lf oc (£l, H) , such that Uk — 1 u. Assume that there exists a dense subset D E ff" 
such that, for any h E D, the sequence {{Puk, h)) is relatively compact in H^" l (n). 
Moreover, let q £ C(0, fC(H)). 

(i) If q — q* and for all (x, £,h) E fl x 5 n_1 x H, one has 
{p(x,£)h = 0) ((q(x)h,h)>0) 
then, for any nonnegative ip E Cq°(0) 

liminf / ip(q(x)uk,Uk)dx > / ip(q(x)u,u)dx. 
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(ii) If for all (x,£, h) G x S n 1 x H, one has 

(p(x,Oh = 0) => ((q(x)h,h)=0) 

then 

(q(x)uk,Uk) converges to (q(x)u,u) inV'(Vl). 

The previous theorem holds also if instead of a differential operator P G OPS m 
we consider the differential operator 

Pu{x) = d a (a a (x)u(x)), 

\a\ <m 

where a G N n , \a\ < m, a a G C(Q, C(H, H^)). 

3. Multidimensional Framework 

We shall restrict our analysis to the case E(W) = to simplify the computations, 
but our results easily extend to the case of E(-) G Lip(M ,R ). 

3.1. Decoupled system. We will consider the following semilinear system of 
equations 

(3.1) W t {x, t) + A(x, D)W(x, t) = B{x, W(x, t)) + D(W(x, t)) 

where t > 0, x G M. d , W G E^. We assume the following hypotheses hold. 
(A.l) B( v ) G C 1 (R d+N ,R N ), £>(•) £ C^K^.R*) 

(A. 2) A(x, 13) G OPS 1 , the system ( |3.1[ ) is hyperbolic, namely the principal symbol 
of A(x, D) is the matrix A(x, £) G AIatxat whose eigenvalues for (x,£) G 
l* 1 x l d , ( / are pure imaginary, 

(A.3) denote by S = span {B(x, W) \ (x, W) G R d+N } then dim S = N - k, < 
k < N. 



Remark 3.1. We point out that (3T) includes the case of the following hyperbolic 
semilinear system 

d 

W t {x, t) + J2 ^j{x)djW(x, t) = B(x, W(x, t)) + D{W(x, t)) 
j=i 

where AAx) € Mnxn, j = 0, . . . , d and for all nonzero vector £ G Mr, the matrix 

j4(x, £) = A;(x) has real eigenvalues. 

i=i 

The aim of this section is to decouple the full system in order to single out the 
conserved quantities from the others. We want to motivate the idea of decoupling 
by considering Carleman's system: 



(3.2) 



/lt + -/l*=4(/2 + /l)(/2-/l) 



f2t fix — —-{fl + - fl)- 

e e- 
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Setting p = fi + and to = we get p, to satisfy 



(3.3) 



pt + m x = 

e 2 TO t + p x = -2pm. 



In this way we decoupled (3.2) in (3.3) and we have isolated p that can be easily 
seen, to be the conserved quantity. Now we want to do the same procedure on 
our system (3.1). Let us consider the hypothesis (A. 3), then there exists a matrix 
P 1 G MkxN such that P I B(W) = and Z 1 = P ! W is the conserved vector. 



Therefore, we can construct an invertible matrix P 



P 1 
p ii 



P" G M 



(N-k)xN- 



Now, for any W G 



wc set 



Z 1 = P T W, 



Z = 

PA{x,C)P- 1 = 



Z 1 
Z u 

P 1 D{P- 1 Z) 

~M xl (x,t) M 12 {x,C) 
M 21 {x,£) M 22 (a;,0 



Z 

QiZ 1 .^ 11 

M(ar,0 M ij (x,D 



11 



P H W, 



P H B(P- 1 Z) 
P II D(P- 1 Z) 
OPM ij (x,£), 



hence by using the previous notations we can rewrite the system (3.1) in the fol- 
lowing decoupled form 



(3.4) 



Z{ + M n (x, D)Z I + M 12 (x, D)Z H = D^Z 1 , Z 11 ) 

Zl 1 + M 21 {x, D)Z I + M 22 {x, D)Z H = Q(x, Z 1 , Z H ) + D II {Z I , Z 11 ), 



where by construction Z 1 = Z J (x, t) GR fe , Z H = Z H (x, t) G R N - k and M n (x, £) G 
M k x k , M 12 (x, G M k x (A r_ fe ) , M 21 (x, i) G M {N ^ k) x fc , M 22 (x, £) G A4 (A r_ fe) x ( , 
D^Z^Z 11 ) G M fe , D II (Z I , Z 11 ) G R^, Q(a;, Z', Z J/ ) G R w " fe . The previous 
transformation does not affect the hyperbolic character of our system. 



3.2. Structural conditions. In order to perform our analysis on system (3J) we 
need the following structural assumption: 

(S.l) M n (x,£) = 0, for any (sc,f) G R d x R d . 

This condition is natural if we consider system ([3.3]). In that case M(x,^) is anti- 
symmetric and is given by 

'0 L 



M(x,0 = 



1 



The same hypothesis is assumed also by Lions and Toscani |27| . It is also implicitly 
contained in the work of Marcati and Rubino |35| . In fact in (3^] the relaxation of 
the following quasilinear nonhomogeneous 2x2 hyperbolic system is considered, 



(3.5) 



( w s + f (w, z) y = 
\z s + g(w,z) y = h{w,z) , 

where y G R, s > 0, with the assumption 

/Ko) = o, 
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which in the linear case is equivalent to (S.l). An analogous condition is assumed 
by Lattanzio and Yong in [p4[ . where they study the singular limits for the initial 
value problem 



\V t + A 3 (W)W Xj {eW)W Xj = , 

3=1 3 = 1 



W(x,0) = W (x-e) 7 

by validating the formal asymptotic approximation in the framework of H s - smooth 
solutions. In 24 they set 



A) 1 A?\ . (Af Af 
A] 1 A h ) A >-\Af Af 



3 3 

and the key structure condition is given by 

A] 1 {0) = 0, for all j. 

In practice the condition (S.l) is essential otherwise the only relaxation process 
would be the trivial one, which relaxes on the null solution. In the case (S.l) is 
not valid, we actually, do not have relaxation from hyperbolic to parabolic systems 
but we have a multiscale phenomena which involves the simultaneous action of dis- 
tinct relaxation mechanisms. To clarify this issue, we follow the formal asymptotic 
approximations of Lattanzio and Yong |]24f , namely we consider the scaled system 

Q(W) 



(3.6) W t + -J2A 1 (x)W, 



e 

3=1 



where W € K, {x,t) G R d x M + , Aj(x) arc smooth N x N matrices with the 
following assumption 

(a) the first k components of Q{W) are zero that is 



Q(W) = 





q(W) 



with q(W) e R N - k . 
Corresponding to this decompostion of Q, we set 



Furthermore we assume: 

(b) q(u, v) = if and only if v = 0, q v (u, 0) is invertible for any u. 
We look for a solution of the form 

oo oo 

O e (x,t) ~^e fc O fe (a;,t) =O {x,t) + J2 sk °k(x,t). 



k=0 fe=l 



Recalling (3.6) we introduce 



d 



R(W) = W t + jJ2Aj(x)W Xj 

e 3=1 



Q(w) 
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Since O e is expected to solve (3.6) we want R(O e ) — 0. Thus it follows 

Q(O ) = o 

d 

J2 A Ax)O 0Xj =Qw(O )O 1 

3=1 

d 

Okt + ^ A j( X )°k+l,Xj = Qw(O )O k+2 
3=1 

The above equations can be rewritten as 

q(u ,v ) = 
d 

£ {A)\x)u 0Xj +Af{x)v 0Xj ) = 
3=1 

d 

} J (Af (x)u 0Xj + A 22 (x)v 0xj ) = q u {uo,Va)ui + q v {u a , v Q )vi 

3=1 



Set O k = 
(3.7) 

(3.8) 
(3.9) 
(3.10) 
(3.11) 



Ukt + ^2 ( j4 ) 1 ( a; )'"fc+i^3 + A ) 2 {x)vk+i, Xj ) = 

3=1 
d 

Ufet (^ 1 (x)M fe+ i, 2 , j . + A 22 ^)^^) = g„(«o,«o)'«fe+2 

3=1 

+ q v (u (h v )v k+ 2- 



From (3.7) and (b) it follows vq = 0, then (3.8) reduces to 

d 

^AfWuta, =0, 

3=1 



so Uoej = 0. By using (b) and (|3.9|) we get = and setting k = in ( 3.10 ) entails 
Mi = 0. Up to now we found ito = vo = U\ = v\ = 0. By assuming inductively 
Up = v p = 0, by using the previous relations and (b) we get, u p+ \ = v p+ ± = 0. 
Hence the formal limit is the null solution. Now, let us consider (3.1), then, by 
using the previous notations and by denoting 



(3.12) 



M(x,D) = OP 



M 12 {x,t) 
M 21 {x,i) M 22 (x,£) 



we can rewrite the system (3.4) in the following form 



(3.13) 



Z t - M(x,D)Z = Q{x,Z\Z 11 )+D{Z 1 ,Z u ) 



where D{Z 1 , Z u ) = (D 1 (Z 1 , Z U ),D II {Z I , Z 11 )). 

We formulate here the hypothesis concerning the existence of a symmetrizer for the 
system ( |3.13j ) in pseudodifferential form (see Taylor p7|), namely 
(A.4) there exists R{x, D) 6 OPS% such that R{x, D)M (x, D>f(R(x, D)M (x, D))*e 
OPSi and its symbol R(x,£) is a positive definite matrix for |£| > 1. 
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The next structure condition regards the existence of a symmetrizer for the system 
( [3.13| ). We assume here a special "block structure" which is natural for strictly 
hyperbolic systems (for instance, in a more complicated framework, see the seminal 
paper of Kreiss JTsj] , Majda and Osher |30| or Ralston The block structure 

follows also for non - strictly hyperbolic systems having constant multiplicity, by a 
general result due to Metivier p?j . Here we are not assuming anyone of the previous 
conditions but directly the "block structure" of the symmetrizer R(x,D). 

(S.2) the symbol of R(x,D) has the following form 



Ru(x,0 o 

Rw(x,C 



where Rit(x,g) G Mkxk, RmfaO G M( N - k )x(N-k) are symmetric positive 
definite matrices and R(x,D)M(x,D) + (R(x, D)M{x, D))* e OPS%. 

Let us remark that in many applications this requirement will be automatically 
satisfyied. 

3.3. Dissipativity condition. In this section we state the assumption on the 
nonhomgeneous term Q(x, Z l , Z 11 ). 

(D) Q(x, Z 1 , Z n ) has the following form 

Q(x, Z 1 , Z H ) = Q (x, Z 1 , Z u ) + Q x {x, Z 1 , Z u ) 
and Q(x, Z ! ,0) = for any (x, Z 1 ) G R d x K fe moreover 

(dl) QofaZ^Z") S C 1 {R N+d ;R N - k ), Q Q „(x, Z 1 , Z n ) is bounded in 

(x, Z 1 , Z 11 ) and [Qq v , R^{x, D)] — 0. There exists Ao > such that for any 

x e R d , (Z 1 , Z u ) e R k x R N - k , Q 0v {x, Z 1 , Z n ) < -X Q I, 

(d2) Q^x.Z^Z 11 ) e C 1 (K Ar+1 ;M Ar - fe ), Q lv (x, Z 1 , Z H ) is bounded in 
(x. Z 1 , Z 11 ) and the operator i?22(s, D)Qiv satisfies 
\\R22{x,D)Q ly \\ c{L 2 ) < Ai, Ai > 0, Ai < A /2. 

Remark 3.2. We splitted the nonhomogeneous term Q(x, Z 1 , Z 11 ) into two terms 
that take into account different dissipativity mechanisms of Q. The former term 
QofcjZ 1 ,Z IT ) cares about the dissipativity of Q, and it commutes with i?22^, D)Qi u . 
The latter Qi(x, Z 1 , Z 11 ) does not commutes but defines with i?22(a;, D) a bounded 
operator dominated by the dissipative part of Qo(x, Z 1 , Z 11 ). 

Remark 3.3. The class of dissipativity terms defined in (D) is not empty. In fact it 
is sufficient to take Q of the following form 

(3.14) Q(x,Z I ,Z II ) = C(x)Z" 

with C{x) G M(N-k)x(N-k) for an y x G R d , [R^ix, D), C(x)} = and there exists 
7 > such that for any x G M d , C(x) < -7/. 
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3.4. Formal analysis of the singular limit. We will analyse the relaxation 
process of the following system 

(3.15) W t + -A(x, D)W = ^B(x, W) + -D(W). 

e e A e 



Following the construction of the previous paragraph we rewrite (3.15) in this way 
(3.16) 

Z{ + -M 12 (x, D)Z n = -DHZ 1 , Z u ) 

£ £ 

Z ? + -M 21 (x, D)Z I + -M 22 (x, D)Z n = \q(x, Z 1 . Z u ) + -D 11 {Z 1 , Z u ). 

By formal asymptotics we are leaded to define 

(3.17) Z\x, t) = U*(x, t), Z n (x, t) = eU n (x, t). 

The previous scaling has an equivalent interpretation as a scaling of the time vari- 



able setting d T = sd t , for more details see |35|. In this way the system ( p.l6[) 
transforms into 

(3.18) 

Ul + M 12 (x,D)U n = -D^U^eU 11 ) 

£ 

£ 2 U t H + M 21 (x, D)U J + sM 22 (x, D)U H = ^Q(x, U 1 , eU H ) + D U (U T , eU H ). 



If we denote by (U 10 , U 110 ) the limit profile as e J, 0, the formal limit system (3.18) 
relaxes to the system 



(3.19) 



Ui° + M 12 (x, D)U IW = Di(U 10 , 0)U IIa 

M 21 (x, D)U m = Q„(s, U 10 , 0)U HO + D H (U 10 , 0), 
where Q u , D l v , denote the derivative respect to the variable irrespectively of Q 



and D 1 . By using the hypothesis (D), the system ( |3.19D transforms into 
U t m + M 12 (x,D)U no = Di(U m ,0)U HO 

U IIQ = Q-\x,U I0 ,Q) [M 21 {x,D)U Ia - D u (U IO ,0)] 
That is equivalent to (by setting U = U 10 ) the second order parabolic system 

(3.20) 

U t + M 12 (x,D)Q-\x,U,0)M 21 (x,D)U = M 12 (x, D)Q- 1 (x 1 U, 0)D n (U, 0) 
+ D I l ,(U,0)Q- 1 (x,U,0) [M 21 (x,D)U - D n {U,0)] . 

In the next sextion we will find sufficient conditions in order to justify rigorously 
this formal analysis. 
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4. A PRIORI ESTIMATES AND CONVERGENCE ANALYSIS 

In this section we consider our rescaled system 
(4.1) 

U{ + M 12 (x,D)U u = -D^U^ell 11 ) 

e 2 U t H + M 21 (x, D)U T + eM 22 (x, D)U H = ^Q(x, U 1 , eU u ) + D 11 ^ 1 , eU n ). 
and we are going to develop the rigorous theory in order to get the relaxed system 



( 3. IS ). We want to show that, as £ J. 0, the solutions of the rescaled system satisfy 

U 1 — ► U 10 strongly in L 2 oc (R d x R+), 

jjil jjiio weakly in L 2 (R d x [0, T]), 

eU H — ► strongly in Lf oc (M. d x R+). 

The basic idea used in this section is the assumption of the existence of a sym- 
metrizer for the system whose symbol has a suitable block structure (condi- 
tions (A. 4) and (S.2)). The pseudodifferential nature of the symmetrizer, as we will 
see later, will allow us to use the properties of such kind of operators (Theorem 2.1, 
Theorem 2.2) in order to obtain "energy" type estimate. 

4.1. A priori estimates. In this section we wish to establish a priori estimates, 
independent of e, for the solution of the system ([O]). To achieve this goal the 
following hypotheses are needed 

(B.l) M 12 (x,0,M 21 (x,0,M 22 (x,0 G S\ 

(B.2) det \(M 21 (x,Z)) T M 21 (x,0~\ ^ 0, 

(B.3) D = (D I {Z I ,Z n ),D II (Z I ,Z 11 )) e C^R^jR* x W N ~ k ), D^Z^O) = 0, 
D H (0, 0) = 0, and Dl is bounded in (Z 1 , Z 11 ), D 11 is a lipschitz function in 
(Z',Z n ), with lipschitz norm a. 



Remark 4.1. Since M 21 {x,^) € M(N~k)xk from elementary linear algebra we de- 

2 



duce condition (B.2) is violated whenever k > ^ 



The next results concerns the "energy" type estimates independent on s, obtained 
via the existence of symmetrizers and via the dissipativity conditions (D). 

Theorem 4.2. Let us consider the solution {U 1 }, {U 11 } of the Cauchy problem for 

,N-k 



the scaled system (gj[). Assume U T (x,0) G [L 2 {R d )] , U H (x,0) G [L 2 (R d )] J 
and the hypotheses (A. 4), (S.2), (B.l), (B.2), (B.3), (D) hold. Then, there exists 
Sq > 0, such that for e G (0, £q), one has 

(i) for all T > 0, there exists M(T) > 0, independent from e, such that 
\\U U \\ L2(M 4 X[0 , T]) < M(T) and sup \\eU n (; i)|| < M{T), 

[0,T] 

(ii) {e 27 !} 1 } is relatively compact in iJ^](R d x R+), 

(iii) {U 1 } is uniformely bounded, with respect to e, in L°° (R + ,L 2 (R d )), namely 

for all T > 0, there exists M(T) > 0, independent from e, such that 

sup || U\;t) || <M(T). 
[0,T] 
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Proof. To simplify our notation, we set 



U 



u 1 



F 



-DUU^eU 11 ) 

e 



Q{x, U 1 , eU 11 ) + D'^U^eU 11 ) 



M e (x,D) = OP 



_£ 

M 12 (x,£) 
,M 21 (x,0 eM 22 {x,i) / 
By using the Theorem 2.2, let R{x, D) be the symmetrizer given by the hypotheses 
(A.4), (S.2), then there exists R = R mod OPS^ 00 , with R > rjl > 0. So we can 
write R= R + T,T G OPS- 00 , where by the block structure of R{x,D), T(x,£) 

(x,£) o 

T22(X,{) 



) . Hence it follows 
: {R{x, D)U e , U £ ) 2 = {R{x, D)U £ t , U e ) 2 + (R(x, D)U e , U e t ) 2 



takes the following form T(x,^) = ( Tii ^q 

d_ 

dV 

= (R(x, D)(M £ {x, D)U + F), U) 2 + {R(x, D)U, M e (x, D)U + F) 2 
= ((R(x, D)M E (x, D) + (R(x, D)M E (x, D))*)U, U) 2 

+ 2(R(x,D)F,U) 2 =Ii+h. 
We estimate separately I\ and I 2 . By applying Theorem 2T and Theorem 2fl and 
by using the hypothesis (A.4) for 1\ we get 

h < 2\\U II \\\\U I \\+e\\U II \\ 2 . 
Let us focus our attention on I 2 , then one has 



(R(x, D)D(U I , eU 11 ), U) 2 + (R(x, D) 







e- 1 Q{x,U I ,eU 



ii 



= (e-'D'iU^eU 11 ), (R n (x, D) + T u (x, D^U 1 )? 



+ {D u (U 1 ,£U 11 ),(R 22 (x,D) +T 22 (x,D))U 11 ) 2 

+ (e^R 22 (x, D)Q{x, U 1 , eU 11 ), U n ) 2 + (e- 1 T 22 (x 1 D)Q(x, U 1 , eU 11 ), U H ) 2 

— h.l + h,2 + -^2,3 + -^2,4- 

To estimate I 2j i and I 2}2 we use hypothesis (B.4) and we get 

l2,i + l2,2<a\\U II \\\\U I \\+ae\\U II \\ 2 . 

Now we turn to 7 2 ,3- Let us denote by Q u = J Q Q v (x, U 1 , e6U II )d6, then we can 
rewrite I 2 ^ in the following way 

h, 3 = {R22{x,D)Qo,U II ,U II ) 2 + (R 22 (x 1 D)Q lu U II ,U II ) 2 = / 2 , 31 +/ 2 , 32 . 

From (dl), it follows 

h,3i = (i4 /2 (^ D)Q 0v U H , Rli 2 (x, D)U h ) 2 = (Q ,i4 /2 (z, D)U u ,Rg\x, D)U H ) 2 

< -Aoll^H 2 . 
By using (d2) we get 



/ 2 ,32 = {R 1 2 , 2 2 {x,D)Q lv U 11 ,R^ 2 {x,D)U 11 ) 2 < X^U^W 2 < 'f\\U u \\ 



Til pl/2/ 



A 



i7||2 



hence 



A 



1/112 



/2,3 = -fll^ i || 
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Now it remains to estimate 1%^. For any 5i > we have 



/ 2 ,4< WT + f \\u u \\ 2 H - 

01 



By adding I\, I2, for all Si > 0, 62 > 0, we get 
d 
It 

c 



^-(R(x, D)U-% W) 2 < (2 + a)||ff"|| \\U X \\ + e(a + 1)\\U H \\ 2 ^||E/"|| 2 



ZZ112 



£1 



< (6 1 +6 2 )\\U n \\ 2 + 



2 + a, 
~2 



||£/ J f + e{a + l)\\U H \\ 2 



If we choose e < 



(a+1) 

(4.2) 

Let us denote by 
(4.3) E{t) 



-yll^f + ^ll^lll-.. 
" and [6\ + S2) < — it follows 



±(R^,D)U s ,Un 2 < ~\\U n \\ 2 + c\\U I f + c\\U II f H ^. 



e 2 \U H (x,t)\ 2 dx+ / |f/ z (a;,f)| da; 



By the second equation of the system (4.1) and by the smoothness of the coefficients 
we get the following estimate 

By integrating (4.2) on [0, f] we obtain the energy E(t) satisfies the following in- 
equality 



(4.4) 25(f) < E(0) ^ 



\U n (x,t)\ 2 dxds + c / E(s)d 



for some c > 0, then by applying the Gronwall's lemma, we have 

(4.5) E(t) < E{0)e ct , 

and 



(4.6) 



JM d 



U H \ 2 dxds < E{0) (t 



We can conclude that for any T > there exists M(T) > 0, indipendent from e, 
such that 

\\U n \\ L 2 {K a x[0tT]) < M(T), sup \\U H (;t)\\ < M(T), sup ||C/'(-,f)|| < M(T). 

[0,T] [0,T] 

In this way we proved (i) and (iii). Let us consider u> relatively compact subset of 



then 



\e 2 U t "\\ H -i M 



sup \(e 2 Uf I ,4>}\= sup 
II^IIhi m =1 ll^llffi^)^ 1 

<e 2 sup (\\U H \\Ut\\)<M(T)e 2 
ll*llffi (t ,)= 1 



e 2 U H (t) t dxdt 



□ 
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Remark 4.3. If there exists a way to control ||Q„||z,j P we do not need any assumption 
(D) but to estimate I 2 ,3 it is sufficient to have for any (x, Z 1 , Z u ) G R d x R k x R N ~ k , 
Q u (x, Z 1 , Z u ) < -XI, A > 0. Indeed 

h, 3 = (R 1 2 / 2 2 (x,D)Q,U II ,R 1 2 / 2 2 (x,D)U II ) 2 

= {Q v Rl! 2 2 {x, D)U n ,Rli 2 (x, D)U H ) 2 + {[R]! 2 (x, D),Q v ]U u , U H ) 2l 



then by taking into account the Proposition |2.3] and the dissipative condition on 
Q v {x, U 1 ,eU n ), we have 

I2fl<-~\\U"\\> + C\\U"\\*H-*- 



which lead to the energy inequality (4.4). 

4.2. Strong convergence analysis. We begin with an immediate consequence of 
(i) and (ii) in the Theorem (p~2|). 

Theorem 4.4. Let us consider the solution {U 11 } of the Cauchy problem 
for system 0>. Assume U 1 (x, 0) G [L 2 (R d )] k , U n (x,0) G [L 2 (R d )] N ~ k and 
moreover the hypotheses (A. 4), (S.2), (B.l), (B.2), (B.3), (D) hold. Then there 
exists U 110 G [L 2 (R d x [0, T])] N ~ k , such that, as £ J. 0, one has (extracting eventually 
a subsequence) 

(4.7) U n ->> U IIQ weakly in L 2 (R d x [0,T]) 

(4.8) eU H — > strongly in L 2 oc (R d x K+) 

(4.9) {e 2 U{'}^0 m H- l (R d xR + ). 



Our next step is to prove the strong convergence for the sequence {U 1 } in the 
norm of L 2 oc (R d x M+). For this porpuse we only use the estimates obtained in the 
previous section. Our main tool in the limit process will be Tartar's and Gerard's 



Theorem 2.4, 1491, 11 



Theorem 4.5. Let us consider the solution {U 1 }, of the Cauchy problem for 
system Assume U*(x,0) G [L 2 {R d )] k , U H (x,0) G [L 2 (R d )] N ~ k and the 

hypotheses (A.4), (S.2), (B.l), (B.2), (B.3),(D) hold. Then there exists U IQ G 
[L 2 (R d x [0,T])] , such that, as e 10, one has (extracting eventually subsequences) 

(4.10) U 1 > U IQ strongly in Lf oc (R d x R+) 



Proof. By using the hypothesis (D), Q^(x,U I ,eU n ) is uniformely bounded in 
L°° then e^Q^x, U 1 ,eU u ) = Q v (x, U 1 , eU")U" is uniformely bounded in L 2 , 
therefore e~ 1 Q v {x, U 1 , elf 11 ) is relatively compact in -ffi~*. In a similar way, 
thanks to the conditions on the function D(J7 / , ell 11 ), we get e~ 1 D I (U 1 , eU 11 ) 
and D 11 (U 1 , eU 11 ) relatively compact in Now by using the identities 

U 1 + M 12 (x,D)U H = ^D^U^eU 11 ) 

M 21 (x, D)U T = -Q(x, U 1 , eU 11 ) + D 11 ^ 1 , eU 11 ) - e 2 U t H - eM 22 (x, D)U H 
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and by taking into account the a priori estimates and the smoothness of the coeffi- 
cients of the system, we can conclude 

f U\ + M 12 (x,D)U IIS 



(4.11) 



M 21 (x, D)U T 



is relatively compact in (H, ~) 



Since the symbols M %J (x,£) are polihomogeneous in the sense of the Section 2, we 
can decompose the operator (x, D) in the following way, 

(4.12) M ij (x, D) = Ml 3 (x, D) + M^(x, D), 

where M[ j (x,D) e OPS 1 , while M$(x,D) 6 OPS . By applying the Theorem 
fo| and by the decomposition ( 4.12 ) we get 

f V{ + Ml 2 (x,D)U u ^ 



(4.13) 



is relatively compact in (H io l) 2 - 



Ml 1 {x,D)U I 

In order to fit into the framework of the Theorem |2.4| we set 
U 1 



P 



u 



11 



dt 



u 1 



Ikxk 



The principal symbol of P is given by 

Ikxk 



for any £ = (£',&) G K d+1 , |£| = 1. We notice that 






Ml 2 (x,D) 




'U 1 ' 


Ml 1 (x,D) 










MP(x : ?) 
M 2 \x,a 







Ml l (x,i')\ = 



for all AeR fc ,/je R N ~ k . 



Now if £' = then ^ ^0 and so A = 0, otherwise if £' ^ by using the hypothesis 
(B.2) we get A = 0. Therefore we get 



(x, £, A, /i) such that p(a;,£) 



A 




A 





(U C {A | A = 0}. 



We take now 



q(x) 



Ikxk 





and for all £ ^ 0, £ = (£o, £') we have 



JV-fe 



= 



implies 



A 




A 


A 


5 


A* 



> = o 



for all Ael l ,/i6 

Now we can apply Theorem (2.4) of Gerard and we conclude that for any 
tp G V(R d x R + ) 



|J7 J (a;, t)\ ip(x,t)dxdt 



U 10 (x,t)\ (f(x,t)dxdt 



where U 10 denotes, in view of Theorem (4.2) the weak limit of U 1 in L 2 (M. d x . 
In this way we obtain 

U 1 



u 



10 



strongly in Lf oc (l 



□ 
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Corollary 4.6. Assume that the hypotheses of Theorems (4-4): (4-') hold, then 
(U 10 ,U 110 ) verifies the following system, in the sense of distributions, 



(4.14) 



Tjio + m 12 (x, D)U H0 = Dl(U 10 , 0)U 



I/O 



M 21 (x,D)U 



io 



,{x 1 U IO ,0)U no + D II (U 10 ,0) 



Proof. By taking into account the regularity hypothesis on Q, the strong conver- 
gence of {U 1 , eU 11 } in Lf oc and the weak convergence of U 11 , we get 



e^Q^x, U^eU 11 ) Q„(s, U 10 , 0)U IIQ weakly in L 
In a similar way we have 



loc ( 



e-'D 1 {U 1 ,eU M ) ^ Dl{U lu ,Q)U 1W weakly in Lf oc (R x R + ). 

There fore we can pass into the limit the other terms and we obtain the relations 
( PI ). □ 



Remark 4.7. We can weakly relax the assumption concerning the polihomogeneity 
of our symbols by assuming directly the decomposition (4.12) 

4.3. Parabolicity condition. Let us restrict our attention to the differential op- 
erator case, namely fll.lD takes the form 

1 d 1 1 

W t + - E A j (x)d j W = — B(x, W) + -D{W). 

In this case the identities ( [I.14| ) become 
d 

Ui° + J2 M}\x)djU no = Di{U m , 0)U 



(4.15) 



no 



2 M 21 (x)d 3 U 10 = Q v {x, U IO ,0)U IIQ + D U (U 10 , 0), 
which is equivalent to write (where we set U = U 10 ) 



(4.16) 



3 = 1 

d 



fc=l 



J2 Mj\x)d j (Q- 1 (x, U, 0)D H (U 1 0)) 

3=1 

d 



+ Di(U,Q)Q- 1 (x,U,0) 



Y d Ml 1 {x)d k U-D II (U,Q) 



Lfc=l 



We want to show, in this simpler case, that (4.16) is parabolic in the sense of 
Petrowski as recalled in the Section 2. Taking into account the notations of Section 
2, one has 

L = M 12 (x,OQu 1 ^>U,0)M 21 (x ) G S 2 
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d d 

where M 12 (x, f ) = i ^ M/ 2 (x)^ and M 21 (a;, £) = i ^ A/* 2 . We are going to 

3=1 3=1 
show the existence of a positive matrix Pq such that Po£ + £*-Po is negative definite. 

By the definition of symmetrizer for the system (3.4) we know the follo wing relations 

between the blocks of R(x, £) and the coefficients of the system (3.4) 

Ru{x,Z)M 12 (x,Z) = (M 21 (x,0) T R22(x,0, 
' '' ' ' i? 22 (x,OM 21 (a;,0 = (M 12 (x,0) T Rn(x,0- 

Now let us denote by 

P = Rn(x,O, 

we have to prove 

P L + L*P <-(3I /3>0. 
By using the relations ( |4.17| ) and the condition (D) it follows 

((P L + L*P )ri,ri) = (R 22 (x,OM 12 (x,OQ- 1 (x,U,0)M 21 (x,^r 1 ,ri) 

+ (AT 21 (x, 0) T (Q^(x, U, 0jf(M 12 (x, 0) T R n (x, OV, V) 
= ((M 21 (x 1 0) T (R22(x,OQ^ 1 (x,U,0)r ] ,7 1 ) 
+ ((Q-^x^.O)) 7 'R 22 (x,Q)M 21 (x,Z)r), V ) 

<-\ \M 21 (x,Ov\ 2 <-(3\vf 
with (3 > 0, this is the notion of parabolicity introduced in the Section 2. 



5. The constant coefficients case 

In this section we want to show how in the case of constant coefficient differential 
semilinear systems our theory can be easily simplified. First we remark that in 
this case we don't need to use pseudodifferential theory because of the constant 
coefficient we can handle them with classical methods. We point out that also in 
this case we will assume the existence of symmetrizers R(x, D) with block structure 
but since the coefficients are constant the principal symbol of R(x, D) depends only 
on the variable £ more exactly it is a homogeneous radial function of degree zero 
of £, R{Cj = R fjlj- j ■ Then the symmetrizers reduce to Fourier multipliers. 

5.1. A priori estimate. Here we consider the following system 



(5.1) 

f d 1 

U 1 + J2 Mj 2 d 3 U H = -D'iU 1 , eU 11 ) 

d d 

e 2 U t H + MfdjU 1 + e £ M 22 d 3 U H = -Q(x, U 1 , eU 11 ) + D 11 ^ 1 , eU 11 ). 

3=1 3=1 £ 

with the same hypotheses (A.4), (S.2), (B.l), (B.2), (B.3), (D) of Section 3.2 and 
Section 4.1, specialized to our simpler framework. Then we have also in this case 
the following theorem 
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Theorem 5.1. Let us consider the solution {U 1 }, {U 11 } of the Cauchy problem for 
the system fiTty. Assume U 1 \x,0) £ [L 2 (R d )] k , U u (x,0) G [L 2 (R d )] N ~ k and that 



hypotheses (A.4), (S.2), (B.l), (B.2), (B.3), (D) hold. Then there exists e > 0, 
such that for s G (0, eo), one has 
(i) for any T > 0, there exists M(T) > 0, independent from e, such that 
\\U n \\ Lm 4 x[0 , T]) < M(T) and sup \\eU u \-,t)\\ < M(T), 

[0,T] 

(ii) {£ 2 C// J } is relatively compact in iJ ; ^](R d x R+), 

(iii) {U 1 } is uniformely bounded, with respect to e, in L°° (M+, L 2 (M d )) , namely 
for any T > 0, there exists M(T) > 0, independent from e, such that 
sup ||tf'(.,t)|| <M(T). 

[O.T] 



Proof. We apply the Fourier transform to the system (5-1) and we multiply in 
L 2 (M. d ) by R(£)U obtaining the following inequality 

+ (R 22 (0U ri ^- 1 QiU I ,eU II )) 2 + (R 22 (OU ri ,LVi(U I ,eU II )) 2 . 



dt 2 



Defining the energy as in ( |4.3| ), taking into account the hypotheses, the properties 
of the symmetrizer and Plancharel theorem it yields the standard energy inequality 



(4.4). The remaining part of the proof follows exactly the same arguments used in 



the previous section so it is omitted. □ 

5.2. Basic ideas on strong convergence. The analysis of strong convergence 
in this case reduces to analyse the convergence of quadratic forms, hence it can be 
obtained via the classical compensated compactness result of Tartar (see E3] , 0] , 
|§ see also §). 

Theorem 5.2. (Tartar's Compensated compactness) 

Let us consider 

1. a bounded open set £1 C M"; 

2. a sequence {1"}%^ V : £1 C R" — ► W n ; 

3. a symmetric matrix : W" 1 — > R m ; 

4. constants a*- fc G R, i = 1, . . . , q, j = 1, . . . , m, k = 1, . . . , n. 
Let us define 

f{a) = (9a, a) , for all a G R m ; 

A= JagM" 1 : 3£gR"\{0}, ^aj^fc 0,i 1.... ,g 
[ i,k 

Assume that 

(a) there exists I G L^ (ft) such that l v — ^ I in L? m (Jl) as v \ oo; 

(b) A l l u = a% jk^~' i = 1, ■ ■ .q are relatively compact in H^ \ (Q); 

(c) /,A = 0; 

(d) there exists / el such that f(l) / m ifte sense o/ measures M(Q). 
Then we have f — f(l). 
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Now we can state our convergence result 

Theorem 5.3. Let {U 1 .U 11 } be the solution of the Cauchy problem for system 
0). Let us assume the hypotheses (A. 4), (S.2),(B.l), (B.2), (B.3), (B.4), (B.6) 
hold, then there exists U 10 <= [L 2 (R d x [0,T])] k and U H0 G [L 2 {R d x [0,T])] N - k , 
such that, as e J, 0, (extracting eventually subsequences) 



(5.2) 
(5.3) 
(5.4) 



U 11 - U H0 



eU 
U 1 



ii 



U 



10 



weakly in L 2 (R d x [0,T]) 
strongly in L 2 loc {R d x K+) 
strongly in L 2 oc (R d x E + ) 



and the limit profile {U 10 , U 110 ) verifies the following system 

d 

U t 10 + M] 2 d 3 U im = Di(U 10 , 0)U J 



(5.5) 



tIIO 



J=l 



M 21 d 3 U m = QviU 10 , 0)U HO + D U (U 10 , 0), 

I 3=1 

in the sense of distribution. 



Proof. ( p.2|) , (|5.3| ) follow from (i) of the Theorem ^jj and the energy estimate 
implies that 



/ d \ 

3=1 

V j=i / 



is relatively compact in (H loc Y 



In order to fit into the framework of Theorem |5.2| we set I s — (U 1 , U 11 ), then the 
characteristic manifold A is given by 



A= {(X,fx) eR k xR N - k | 3£ e R d+1 \ {0} , £?(£, A) = 0} 



where 



\ 



e A+ [Y^MfiAn 



vJ=l 



v v=* 



By using the hypothesis (B.4) and by defining O = ^ ( 7 k ■* k we have /(A) 
A T 0A and, of course, /u = 0. We apply the Theorem 5.2 to show 



70\2 



{U 1 ) 2 - (U w ) 



and finally 



U 1 



u 



10 



in the sense of measure 



strongly in L 2 oc (R d x 



So we can pass to the limit into the nonlinear terms of the system ( |5.lf ). 



□ 
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Remark 5.4. An equivalent form of system ( p.5| ) is given by (setting U m = U) 

d / d \ 



(5.6) 



Ut+Y, M ? d i Qu 1 ^ U, 0) Y M^dkU 



3=1 



fc=l 



X] M T d o (Qu 1 ( x > U : 0)D u (U, 0)) 

3=1 

d 



+ Dl(U I ,0)Q-\x,U,0) 



Y M k ld kU - D n (U,0) 



The proof of the parabolicity of (5.6) follows the same arguments used in the vari- 
able coefficients case. 



6. Approximation of given parabolic systems 

In this last section we want to apply the theory of the previous ones to approx- 
imate a generic given parabolic system, provided the Petrowski condition hold. In 
fact here we reconstruct a parabolic system by means of a suitable larger semilinear 
hyperbolic system that relaxes on it. Two very large classes of parabolic systems 
will be taken in consideration. The former are quasilinear parabolic systems in 
divergence form and the latter are the so called "Reaction- diffusion" systems. 

6.1. Quasilinear case. We consider now the following quasilinear system in di- 
vergence form 



(6.1) 



3 = 1 



where x G M. d , t G R+, U — U(x, t) G R k , N = (d+l)k ( then k < f ). Let us denote 
by F(Z7) G M kxd , (F(U))i = F^U) for any % = 1,... ,d and B(U) G M kdxkd , 
(M(U))ij = Bij(U) for any i, j — 1, . . . ,d, then we assume. 

(C.l) By(-) eC 1 {R k ;M kxd ),for;myi,j = l,... , d and B tJ (U)^ > c \\\ 2 I 
for any A G R d (strong parabolicity), moreover B _1 (i7) is bounded on U, 

(C.2) F. L (-) G C 1 (R k ;M kxd ), for any i = 1,... ,d, F(0) = and B- 1 ([/)F(C7) is 
lipschitz on U. 

(C.3) G(-) G Lip(R k ;R k ). 

We have the following theorem 



Theorem 6.1. Let us consider the system (6.1), let us suppose hypotheses (C.l), 
(C.2), (C.3) holds, then the solutions of the system 



1 



fii 



(6.2) 



Z\ + - div Z 1 

+ ^DZ 1 = —B-^Z^Z" + ^M- 1 {Z I )¥(Z I ), 



where (x, t) G K d x M+, Z 1 = Z I {x,t) G M fe , Z 7/ = Z n (x, t) G M kxd , approximate 
the system (6.1) in the sense of Theorem 4^4 an d the Theorem J^.L . 
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Proof. It can be easily shown that ( |6.2| ) is an hyperbolic symmetric system. By 
rescaling the variables, as in (3.17), the system (|6.2|) transforms into 



(6.3) 



U 1 + divV 11 = GiU 1 ) 



e 2 W 



DU 1 



ii 



^{U^WiU 1 ). 



The conditions (B.l), (B.2) of the Section 4 are satisfied by the system fl3.3p.Now 
by setting D II {Z I , Z H ) = M- 1 {Z I )¥(Z I ) we can easily verify (B.3) and by denoting 
Q(x, Z 1 ,7L U ) = -B- 1 (Z 7 )Z ZJ we get obviously the condition (D). We can apply 
the Theorem 4J3 and we obtain that the solutions to ( |6.3| ) satisfy, as e j (setting 
U 10 = U) 

U t + div(¥{U) + M(U)DU) = G{U). 

□ 

6.2. Reaction - Diffusion type systems. Let us consider the following reaction- 
diffusion type system 



(6.4) 



u t = A jM (x)d 3 d k u + f(u) 

j,k=l 



where x G R d , t G K+, U = U(x,t) £ K fe . We make the following hypotheses. 

(D.l) Aj, fc (.) G C°°{R d ;Mkxk), for any j,k = l,...d and £)J jfc=1 Aj^XjXk > 

co|A| 2 /, for any A G M. d (strong parabolicity) . 
(D.2) /(•) G Lip(R k ;R k ). 

In order to approximate system (5.4), we define the following linear operator; 

Bj(x) : R kd — >R kd , 

d 

B 3 (x)(W) = A jtk (x)W k for any j = 1, . . . , d 

fc=i 



Now system (6.4) can be written in the equivalent form 



U t =Y J B 1 {d J DU) + f{U). 
i=i 

We have the following theorem. 

Theorem 6.2. Let us consider the system §6.4\ ), let us suppose the hypotheses 
(D.l), (D.2) hold, then the solutions of the system 



(6.5) 



Z' + ^B^Z 11 = f(Z I ) 

i=i 
d 

Zi 1 + - ^ BfdjZ 1 = — 5 A{x)Z 



ii 



where (x,t) G R d x R + , Z 1 = Z J (a:,t) G 



7ll 



Z n (x,t) G R fcd , A(» G 



M kdxkd, (A (x))j,k — Aj,k(x) approximate the system \6.4\) ™ ^ e sense of the 
Theorem 4-4 an d the Theorem \ 
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Proof. The system (6.5) is symmetric and hyperbolic. By rescaling the variables as 
in (3.17) it transforms into 



(6.6) 



Ui + '£B j d j lJ 11 = fill 1 ) 



j'=i 



J2 B J d 3 U ' = -A(a?)U- 



ii 



By setting Q{x, Z 1 ,'L 11 ) = —A(x)'Z 11 , the dissipativity conditions (D) together 
with the hypotheses (B.l) and (B.2) of Section 4 are immediately satisfied. Now, 



we can apply the Theorems [4.4| and 4^5 , then since ~^^Bj djU 1 = A(x)DU I , as 

3=1 

□ 



e I 0, the solution of (6.6) satisfy the system (6.4) 



Remark 6.3. The Theorem 6.2 can be applied, with slight modifications, to the 
more general case of Petrowski parabolic systems in the sense recalled in the Section 
2. In fact let us consider the following system 



(6.7) 



d 

u t =J2 c^ k {x)d 3 d k u + G(u) 



where x £ M. d , t £ R + , U = U(x,t) £ M. k and the matrices Cj_k satisfy the Petrowski 
parabolicity condition given in the Section 2. Let us consider the matrix Pq from 
the previous definition, denote by 

W = P U 

then, the system ( |6.7D transforms into a strongly parabolic system, so we can apply 
the previous Theorem 3.2. 



To approximate the system (6.4) we can also follow a different approach by using 
the pseudodifferential theory. In fact, let us denote by 



(6.8) 



(el 11 , the principal symbol of ([O]) and by A(x, D) = OPA(x, £). Now we set 

(6.9) B(x,D)=OP[(-A(x,0) 1/2 } 
We have the following theorem. 



Theorem 6.4. Let us consider the system ^6.^ ), suppose that hypotheses (D.l), 
(D.2) hold, then the solutions of the system 



(6.10) 



Z( + -B(x,D)Z H = f(Z J ) 
Zl 1 --B(x,D)Z I = -\z IJ , 



where (x,t) £ R d x R + , Z 1 = Z I {x,t) £ M fe , Z n = Z n {x,t) £ R k , N = 2k 



approximate the system (6.4 ) in the sense of the Theorem {(-4 and the Theorem 
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Proof. By rescaling the variables as in ( |3.17 ) the system (3.4) transforms into 

U{ + B(x,D)U u = fiU 1 ) 



(6.11) 



e 2 U t n -BfaD)!/ 1 = -U 



ii 



By using (D.l) we have ( 6.11 ) is an hyperbolic system. Let us denote by 
M 12 (x,D) = B(x,D), M 21 (x,D) = -B(x,D), by using the hypothesis(C.l), the 
conditions (B.l) and (B.2) are satisfied. It can also be easily verified that the 
symmetrizer of ( 6.1 0| ) is given by the matrix Inxn- Finally we have to verify the 



dissipativity condition (D). Setting Q(x, Z , Z ) = —Z we get the hypothesis 



(D) is satisfied. Now we can apply Theorem [4.4| and Theorem 4.5 and we obtain 
that the solution of (6.11) satisfy as e \ (by setting U 10 = U) 



U t -A{x,D)U = 0. 

□ 
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